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ABSTRACT 

(N ■ 

^ Context. A variety of formation scenarios have been proposed to explain the diversity of properties observed in bulges. Studying their intrinsic 
shape can help to constraint the dominant mechanisms at the epochs of their assembly. 

Aims. The structural parameters of a magnitude-limited sample of 148 unbarred SO-Sb galaxies were derived in order to study the correlations 
\l ' between bulges and disks, as well as the probability distribution function of the intrinsic equatorial eiiipticity of bulges. 

'^"J Methods. We present a new fitting algorithm (GASP2D) to perform two-dimensional photometric decomposition of the galaxy surface- 
brightness distribution. This was assumed to be the sum of the contribution of a bulge and disk component characterized by elliptical and 

T-H concentric isophotes with constant (but possibly dilferent) eiiipticity and position angles. Bulge and disk parameters of the sample galaxies 
were derived from the y-band images, which were available in the Two Micron All Sky Survey. The probability distribution function of the 
equatorial eiiipticity of the bulges was derived from the distribution of the observed ellipticities of bulges and misalignments between bulges 
^ and disks. 

• ^ Results. Strong correlations between the bulge and disk parameters were found. About 80% of bulges in unbarred lenticular and early-to- 
intermediate spiral galaxies are not oblate but triaxial ellipsoids. Their mean axial ratio in the equatorial plane is (B/A) = 0.85. Their probabil- 
^ , ity distribution function is not significantly dependent on morphology, light concentration or luminosity. The possible presence of nuclear bars 
■ - - I does not influence our results. 

Conclusions. The interplay between bulge and disk parameters favors scenarios in which bulges have assembled from mergers and/or have 
grown over long times through disk secular evolution. However, all these mechanisms have to be tested against the derived distribution of bulge 
intrinsic ellipticities. 

Key words, galaxies: bulges - galaxies: formation - galaxies: fundamental parameters - galaxies: photometry - galaxies: statistics - galaxies: 
structure 



1. Introduction 

The relative prominence of galactic bulges with respect to their 
disks is important in the definition of galaxy types. Therefore, 
understanding the formation of bulges is key to understanding 
the origin of the Hubble sequence. 

Bulges are di verse and h ete rogeneous (see 



the reviews by iKorinendv 19931 : IWvse et al. 1997 



Kormendv & Kennicutt 20041) . The big bulges of lenticu- 
lars and early-type spirals are similar to low-luminosity 
elliptical galaxies. Their surface-brightness radial profiles gen- 
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erally follow a De Va ucouleurs law (lA ndredak is et al. 1995 



Mollenhoff &Hei dt 2001, hereafter 
appear rounder 



ICaroUo et al. 1998 

MHOl). The majority of these bulg es 
than their associated disks (iKent 198 5V Their kinematical 
properties are well described by dynamical models of rota- 
tionally flattened oblate spheroi ds with little or no anisotropy 



dKormendy & Iflingwort h 1982 : 
[Cappellari et al. 2006). They have photometrical 



Dayies & lUingworfli 1983 



matical properties, 
(FP) correlation 



which satisfy the 



Burstein et al. 1997 



Bender et al. 199: 
Aguerri et al. 2005 



and kine- 
fundamental plane 



Bender etal. 1993 



On the contrary, 
the small bulges of late-type spiral galaxies seems to be rem- 
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iniscent of disks. Their surface-brigh tness radial profiles have 
an almost exponential fall o ff ( An dredakis & Sanders 1994 ; 
de Jong 19961: iMacArthur et al. 2003h . In some cases they 



have apparent flatteni ngs that are similar or even larger than 
their asso ciated disks ( Fathi & Peletier 2003b and r otate as fast 
as disks ( Kormendy 1993 : iKormendv et al. 2002 ). Late-type 
bulges deviate from the FP (Carollo 1999). 

Different formation mechanisms (or at least a variety of 
dominant mechanisms at the epochs of star formation and 
mass assembly) were proposed to explain the variety of prop- 
erties observed in bulges. Some of these formation processes 
are rapid. They include early formation from the dissipa- 
tive collapse of protogalactic gas clo uds (Eggenetal. 1962; 
Sandage 199dHGilmore & Wvse 1998HMerlin & Chiosi 2006) 



or later assem bly fro m mergers between pre-existi ng disks 
(iKauffmann 1 996; Bau gh et al. 1996l;ICole et al. 2000l) . In both 
scenarios the disk forms after the bulge as a consequence of ei- 
ther a long star-formation time compared to the collapse time 
or a re-accretion around the newly formed bulge. 

Bulges can also grow over long timescales 
through the disk secular evolution driven by bars 
and/or environmental effects. Bars are present in 
more than half of disk galaxies in the local universe 
(lEskridge et al. 2000l: Menendez- Delmestre et al. 2 00?') 

and out to z ~ 1 ( Eknegr een et al. 2004r Jogee et al. 2 004). 
They are efficient mechanisms for driving gas inward to 
the galacti c center and feed th e galactic supermassive black 
hole (see Corsini et al. 2003a[ and references therein). In 
addition , bar dissolution due to the growth of a central 
mass ( Pfenniger & Norman 1990|). scatte ring of disk stars 
at vertical resonances (ICombes et al. 19 90). and coher- 



ent bending of t he bar perpendicular to the disk plane 
dRaha et al. 1991^, 'Pebattista et al. 2004 lAthanassoula 20051; 
[Martinez- Valpuesta et al. 2006) are efficient mechanisms 
in building central bulge-like structures, the so-called 
boxy/peanut bulges. Moreover, the growth of the bulge out 
of disk material may also be externally triggered by satellite 
accretion during m i nor merging events (Searle & Zinn 1978 : 



lAgue rri et al. 2001'; 'E liche-Moral et al. 2006^ and gas infall 
dThakar & Ryden 19981^ 

Traditionally, the study of the relations between the struc- 
tural parameters of the galaxies have been used to understand 
the bulge formation processes, e.g., the correlation between the 
bulge effective radius and the scale length of the disk in many 
galaxy samples has always been interpreted as an indication 
that bulges were forme d by secular evolution of their disks (see 



MacArthur et al. 2003b . However, one piece lost in this study 



is the three-dimensional shape of the bulges. By studying this, 
one might be able to provide the relative importance of rapid 
and slow processes in assembling the dense central components 
of disk galaxies. A statistical study can provide a crucial piece 
of information for testing the results of numerical simulations 
of bulge formation for different galaxy type along the morpho- 
logical sequence. 

In this paper, we analyze a sample of unbaiTed early- 
type disk galaxies to derive the intrinsic ellipticity of their 
bulges in the galactic plane. The twisting of bulge isophotes 
dLindblad 1956t IZaritskv & Lo 198^) and misalignment be- 



tween the major axes of the bulge and disk (iBertola et al. 199 lb 
are not possible if the bulge and disk are both oblate. 
Therefore, they were interpreted as a signature of bulge triax- 
iality. This idea is sup ported by the presence of non-circular 
gas motions (e.g.. iGerhard & Vietri J98l iBertola et al. 19891; 
Gerhard et al. 1989t lBerman200lb and a velocity gradi- 
ent along the galaxy minor axis ( Corsini et al. 2003bl ; 



Coccato et al. 2004t ICoccato et al. 20051) . We improve the pre- 
vious works in several aspects. First, we use near-infrared 
images to map the distribution of the mass-carrying evolved 
stars and avoid contamination of dust and bright young stars. 
Second, we retrieve the structural parameters of the bulge and 
disk by applying a new algorithm for two-dimensional photo- 
metric decomposition of the observed surface-brightness distri- 
bution. Finally, we obtain the probability distribution function 
(PDF) of the intrinsic equatorial ellipticity of bulges by using 
a new mathematical treatment of the equations describing their 
three-dimensional shape. 

The paper is organized as follow. The selection criteria of 
our sample galaxies and the analysis of their near-infrared im- 
ages are described in Sect.|2l Our new photometric decomposi- 
tion method for deriving the structural parameters of the bulge 
and disk by analyzing the two-dimensional surface brightness 
distribution of galaxies is presented in Sect. [3] The correlations 
between the structural parameters of the sample galaxies are 
discussed in Sect. ID The PDF of intrinsic equatorial ellipticity 
of the studied bulges is derived in Sect.|5] Our conclusions and 
a summary of the results are given in Sect.|6] 

2. Sample selection and data acquisition 

Our objective was to select a well-defined complete sam- 
ple of nearby unbaiTed disk galaxies to study in a sys- 
tematic way the photometric properties of their structural 
components. Since these properties are strongly dependent 
on the dominating stellar population at the observed wave- 
length, it is preferable to consider near-infrared images to 
map the mass-carrying evolved stars and avoid contami- 
nation due to dust and bright young stars (e.g., MHOl). 
The complete sa mple is drawn fro m the Extended Source 
Catalogue (XSC) djarrett et al. 2000b of the Two Micron All 
Sky Survey (2MASS) dSkrutskie et al. 20061) . Our sample con- 
sists of galaxies that meet the following requirements: (1) 
Hubble type classification from SO to Sb (-3 < HT < 3; de 
Vaucouleurs et al. 1991, hereafter RC3) to ensure that bulges 
are fully resolved in 2MASS images; (2) unbarred classifi- 
cation in RC3; (3) total 7-band magnitude Jj < 10 mag 
(2MASS/XSC); (4) inchnation / < 65° (RC3) to measure the 
misalignment between the position angle of the bulge and disk; 
(5) Galactic latitude \ be \> 30° (RC3) to minimize both 
Galactic extinction and contamination due to Galactic fore- 
ground stars. 

We ended up with a sample of 184 bona-fide unbarred 
galaxies. We retrieved their 2MASS 7-band images from the 
NASA/IPAC Infrared Science Archive. The galaxy images 
were reduced and flux calibrated with th e standard 2MAS S 
extended source processor GALWORKS dJarrett et al. 20001) . 
Images have a typical field of view of few arc-minutes and a 



J. Mendez-Abreu et al.: Structural properties of disk galaxies 



3 



-3 -2 -1 



llllllllnn.n 

2000 4000 6000 8000 
V3K (km s-i) 



Fig. 1. Distribution of the sample galaxies over the Hubble 
types (left panel) and radial velocities with respect to the 3K 
background (right panel). 

spatial scale of 1" pixel '. They were obtained with an aver- 
age seeing FWHM~3'.'l as measured by fitting a circular two- 
dimensional Gaussian to the field stars. 

After a visual inspection of the images, we realized that 
some of the sample galaxies were not suitable for our study. We 
rejected paired and interacting objects as well as those galaxies 
that resulted in being barred after performing the photometric 
decomposition (see Sect. |3]l. Therefore, the final sample pre- 
sented in this paper contains 148 galaxies (90 lenticular and 58 
early-type spiral galaxies). A compilation of their main proper- 
ties is given in Table [3] Figure [T] shows the distribution of the 
sample galaxies over the Hubble types. The lenticular galaxies 
are predominant over the spirals due to our magnitude selec- 
tion, which favors red galaxies. Moreover, we show the distri- 
bution of radial velocities of the sample galaxies with respect 
to the 3K background. The mean radial velocity is 2000 km s"' 
(corresponding to a distance of 27 Mpc by assuming Ho = 75 
km s"' Mpc"'), but we include galaxies as far as 8500 km s"' 
(113 Mpc) because the sample is magnitude limited. 

Tonry et al. (2001) derived the distance of 30 galaxies of 
our sample from the measurement of their surface brightness 
fluctuations. The difference between the distances obtained 
from radial velocities and those derived from surface bright- 
ness fluctuations was calculated for all these galaxies. The stan- 
dard deviation of the distance differences is 5 Mpc and it was 
assumed as being a typical distance error. For the 4 common 
galaxies in the Virgo cluster it is 2 Mpc. 

3. Two-dimensional bulge-disit parametric 
decomposition 

Conventional bulge-disk decompositions based on elliptically 
averaged surface-brightness pro files usually do no t take into 
account the intrins ic shapes (e.g., Prieto et al. 200 ll) or the po- 
sition angle (e.g.. iTrujillo et al. 2001ch of the bulge and disk 
com ponents, which can produ ce systematic errors in the results 
(e.g.. lBvun & Freeman 1995h . 



For this reason a number of two-dimensional paramet- 
ric decomposition techniques have been developed in the last 
years. As examples we may point out the algorithms devel- 
oped by Simard (1998, GIM2D), Peng et al. (2002, GALFIT), 
de Souza et al. (2004, BUDDA) and Pignatelli et al. (2006, 
GASPHOT). These methods were developed to solve diff'er- 
ent problems of galaxy decomposition when fitting the two- 
dimensional galaxy surface-brightness distribution. They use 
diff'erent functions to parametrize the galaxy component and 
diff'erent minimizations routines to perform the fit. 

In this paper we present our new decomposition algorithm 
named GASP2D (GAlaxy Surface Photometry 2 Dimensional 
Decomposition). The code works like GIM2D and GASPHOT 
in minimizing the interaction with the user. It works in an au- 
tomatical way to be more efficient when dealing with a large 
amount of galaxies. However, like GALFIT and BUDDA it 
also adopts a Leverberg-Marquard algorithm to fit the two- 
dimensional surface-brightness distribution of the galaxy. This 
reduces the amount of computational time needed to obtain a 
robust and reliable estimate of the galaxy structural parameters. 

In the present work, we show the first version of the code. 
We assume that the galaxy can be modeled with only two com- 
ponents, the bulge and the disk. In a forthcoming paper, we will 
show an improved version of GASP2D with the possibility to 
fit other galaxy components, like bars. 

3.1. Photometric model 

We assumed the galaxy surface-brightness distribution to be the 
sum of the contribution of a bulge and disk component. Both 
of them are characterized by elliptical and concentric isophotes 
with constant (but possibly different) ellipticity and position 
angles. 

Let (^, J], be the Cartesian coordinates with the origin in 
the galaxy center, the ^-axis parallel to the direction of right 
ascension and pointing westward, the 77- axis parallel to the di- 
rection of declination and pointing northward, and the f-axis 
along the line-of-sight and pointing toward the observer The 
plane of the sky is confined to the (^, ij) plane. 

We adopted the Sersic law (ISersic 1968b to describe 
the surface brightness of the bulge component. The Sersic 
law has been extensively used in the literature for mod- 
eling the surface brightness profiles of galaxies. For in- 
stance, it has been use d to model the surface bright- 



ness of elliptical galaxies (IGraham & Guzman 2 003). bulges 
of early and l a te type galaxies (lAndredakis et al. 19951 
Prieto etal. 2OO1I lAguerri et al. 2004 Jvlollenhoff' 200i), the 
low surface brightne s s host galaxy of bl ue compact galax- 
ies ( Caon et al. 2005 ;_ Amorm et al. 20071) and dwarf ellip- 
tical galaxies (iBinggeli & Jerjen 1998t lAguerri et al. 2005bt 
Graham & Guzman 2003h . It is given by 



4(^,7,) = /elO-^»[(^/^=>""-'], 



(1) 



where re, /g, and n are the effective (or half-light) radius, the 
surface brightness at r^, and a shape parameter describing the 
curvature of the surface-brightness profile, respectively. The 
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value of b„ is coupled to n so that half of the total lumi- 
nosity of the bulge is within and can be approximated as 
bn = 0.868 n - 0.142 dCaon et al. 1993h . Bulge isophotes are 
ellipses centered on (^o, ^o) with constant position angle PAb 
and constant ellipticity e\, = I - q^. The radius is given by 



rb = [(-(^-^o)sinPAb + (77- 77o)cosPAb)^+ 
- ((f - fo) cos PAb + (77 - 770) sin PAb)V?b] 



1/2 



(2) 



We adopted the exponential law dFreeman 19701) to de- 
scribe the surface brightness of the disk component 



-1-6 /h 



(3) 



where /q and h are the central surface brightness and scale- 
length of the disk, respectively. Disk isophotes are ellipses cen- 
tered on (^0, 770) with constant position angle PAd and constant 
ellipticity - I - qd- Disk inclination is / - arccos (^d)- The 
radius is given by 



rd = [(-(^-^o)sinPAd + (77 - 77o)cosPAd)^-H 
- ((^ - ^0) cos PAd + (77 - 770) sin PAif/ql] 



1/2 



(4) 



To derive the coordinates (^0, 770) of the galaxy center and 
the photometric parameters of the bulge (/g, r^, n, PAb, and ^b) 
and disk (/q, h, PAd, and qi) we fitted iteratively a model of 
the surface brightness /m(^, ifj-hi^, ifj+hi^, 77) to the obser- 
vations using a non-linear least-squares minimization method. 
It was based on t he robust Levenberg-Marquardt method (e.g.. 
Press et al. 1996h implemented by More et al. (1980). The ac- 
tual computation was done using the MPFIT algorithm im- 
plemented by C. B. Markwardt under the IDL environment. 
MPFIT allows the user to keep constant or impose boundary 
constraints on any parameter during the fitting process. 

For each pixel (^, 77), the observed galaxy photon counts 
/g(f, 77) are compared with those predicted from the model 
/m(f, 77). Each pixel is weighted according to the variance of 
its total observed photon counts due to the contribution of both 
galaxy and sky, and determined assuming photon noise limita- 
tion by taking into account the detector readout noise (RON). 
Therefore, the to be minimized can be written as 



X 



[/„(^, 77)-/g(^, 77)]^ 
^^^g(^>'7) + /s(^,77) + RON 



^_ 1 1 ^ s 



2 ' 



(5) 



with ^ and 77 ranging over the whole N x M pixel image. 

An important point to consider here is the weight function 
used to calculate thex^'^ some authors claim that is better to as- 
sign to each pixel an u ncertainty given by the Poissonian noise 
(e.g.. |Peng et al. 2002h while others adopt constant weights to 
obtain better results (e.g., MHOl; de Souza et al. 2004). Both 
possibilities were implemented in the fitting algorithm. We 
adopted the poissonian weights after extensive testing with ar- 
tificial galaxies. 



' The updated version of this code 
|http://cow.physics.wisc.edu~craigm/idl/idl.html 



is available on 



3.2. Seeing effects 

The ground-based images are affected by seeing, which scat- 
ters the light of the objects and produces a loss of spatial res- 
olution. This is particularly critical in the central regions of 
galaxies, where the slope of the radial surface brightness pro- 
file is steeper. Since the bulge contribution dominates the sur- 
face brightness distribution at small radii, seeing mostly affects 
bulge structural parameters. Seeing effects on the scale parame- 
ters of a Sersic surface brightness profile have been extensively 
discussed by Trujillo et al. (2001a,b). 

During each iteration of the fitting algorithm, the seeing 
effects were taken into account by convolving the model im- 
age with a circular two-dimensional Gaussian point spread 
function (PSF). The PSF FWHM was chosen to match the 
one measured from the foreground stars in the field of the 
2MASS galaxy image. The convolution was performed in the 
Fourie r domain using t he fast Fourier transform (FFT) algo- 
rithm dPress et al. 19961) before thex^ calculation. Our code al- 
lows us to introduce also a Moffat or a star image to reproduce 
the PSF. We tested that the results are not improved by adopt- 
ing a circular two-dimensional Moffat PSF or by computing the 
convolution integrals. 

3.3. Technical procedure of the fit 

Since the fitting algorithm is based on a x^ minimization, it 
is important to adopt initial trials for free parameters as close 
as possible to their actual values. This would ensure that the 
iteration procedure does not just stop on a local minimum of 
the x^ distribution. To this aim we proceeded through different 
steps. 

First, the photo metric package SExtractor (see 



iBertin & Arnouts 19961 for details) was used to measure 
position, magnitude and ellipticity of the sources (e.g., fore- 
ground stars, background and companion galaxies, as well as 
bad pixels) in the images. 

We then derived the elliptically-averaged radial profiles 
of the surface brightness, ellipticity and position angle of the 
galaxy. We fitted ellipses to the galaxy isophotes with the 
ELLIPSE task in After masking the spurious sources 

using the parameters provided by SExtractor, we fitted ellipses 
centered on the position of the galaxy center. This could be 
estimated by either a visual inspection of the image or with 
SExtractor. The coordinates of the galaxy center were adopted 
as initial trials for (^0, 770) in the two-dimensional fit. 

In the third step we derived some of the trial values 
by performing a standard one dimensional decomposition 
technique similar to that adopted by several authors (e.g., 
iKormendy 1977 ; Prieto et al. 2001 ). We began by fitting an ex- 
ponential law to the radial surface-brightness profile at large 
radii, where the light distribution of the galaxy is expected 
to be dominated by the disk contribution. The central surface 
brightness and scalelength of the fitted exponential profile were 
adopted as initial trials for /q and h, respectively. The fitted pro- 



- IRAF is distributed by NOAO, which is operated by AURA Inc., 
under contract with the National Science Foundation. 
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file was extrapolated at small radii and then subtracted from the 
observed radial surface-brightness profile. The residual radial 
surface-brightness profile was assumed to be a first estimate of 
the light distribution of the bulge. We fitted it with a Sersic law 
by assuming the bulge shape parameter to be 0.5, 1, 1.5, ...,6. 
The bulge effective radius, effective surface brightness, and 
shape parameter that (together with the disk parameters) gave 
the best fit to the radial surface-brightness profile were adopted 
as initial trials for re, /e, and n, respectively. 

We also obtained the initial trials for ellipticity and posi- 
tion angle of the disk and bulge, respectively. The values for 
and PAd were found by averaging the values in the outermost 
portion of the radial profiles of ellipticity and position angle. 
The initial trials for q\, and PAb were estimated by interpolat- 
ing at Te the radial profiles of the ellipticity and position angle, 
respectively. 

Finally, the initial guesses were adopted to initialize the 
non-linear least-squares fit to galaxy image, where all the pa- 
rameters, including n, were allowed to vary. A convergent 
model was reached when the had a minimum and the rela- 
tive change of the;i'^ between the iterations was less than 10"^. 
A model of the galaxy surface-brightness distribution was built 
using the fitted parameters. It was convolved with the adopted 
circular two-dimensional Gaussian PSF and subtracted from 
the observed image to obtain a residual image. In order to en- 
sure the minimum in the ;^'^-space found in this first iteration, 
we perform two more iterations. In these iterations, all the pix- 
els and/or regions of the residual image with values greater or 
less than a fixed threshold, controlled by the user, were rejected. 
Those regions were masked out and the fit was repeated assum- 
ing, as initial trials for the free parameters, the values obtained 
in the previous iteration. These kind of masks are usually useful 
when galaxies have other prominent structures, different from 
the bulge and disk (e.g. spiral arms and dust lanes), which can 
affect the fitted parameters, solving the problem in an auto- 
matic way. We found that after three iterations the algorithm 
converges and the parameters do not change. 

3.4. Test on simulated galaxies 

To test the reliability and accuracy of our two-dimensional 
technique for bulge-disk decomposition, we carried out exten- 
sive simulations on a large set of artificial disk galaxies. We 
generated 1000 images of galaxies with a Sersic bulge and an 
exponential disk. The central surface-brightness, scalelength, 
and apparent axial ratios of the bulge and disk of the artificial 
galaxies were randomly chosen in the range of values observed 
in the /-band by MHOl for a sample of 40 bright spiral galax- 
ies. The adopted ranges were 

1 < re < 3 kpc 0.4 < < 0.9 0.5 < n < 6 (6) 
for the bulge parameters, and 

2 < /! < 5 kpc 0.4 < < 0.9 (7) 

for the disk parameters. The parameters of the artificial galaxies 
also have to satisfy the following conditions 

qd<qb < B/T < 0.8 S < Jt < 12 mag. (8) 



All the simulated galaxies were assumed to be at a distance 
of 30 Mpc. This is the average distance of the galaxies of our 
sample and it corresponds to a scale of 145 pc arcsec"^ The 
pixel scale used was 1 arcsec pixel \ and the CCD gain and 
RON were 8 e" ADU ' and 40 e , respectively, in order to 
mimic the instrumental setup of the 2MASS data. Finally, a 
background level and photon noise were added to the result- 
ing images to yield a signal-to-noise ratio similar to that of the 
available 2MASS images. 

The two-dimensional parametric decomposition was ap- 
plied to analyze the images of the artificial galaxies as if they 
were real. Errors on the fitted parameters were estimated by 
comparing the input pi and output po values. Relative errors 
{I- Pi I Po) were assumed to be normally distributed, with mean 
and standard deviation corresponding to the systematic and 
typical error on the relevant parameter, respectively. 

The results of the simulations are given in Table [1] In Run 
1 we built the artificial galaxies by assuming the correct values 
of PSF FWHM and sky level, so only errors due to the Poisson 
noise are studied. The mean relative errors on the fitted param- 
eters are smaller than 0.01 (absolute value) and their standard 
deviations are smaller than 0.02 (absolute value) for all galaxies 
with < 10 mag, proving the reliability of our derived struc- 
tural parameters. Relative errors increase for fainter galaxies, 
which are not included in our sample. 

Systematic errors given by a wrong estimation of PSF and 
sky level are the most significant contributors to the error bud- 
get. To understand how a typical error in the measurement of 
the PSF FWHM affects our results, we analyzed the artificial 
galaxies by adopting the correct sky level and a PSF FWHM 
that was 2% larger (Run 2) or smaller (Run 3) than the ac- 
tual one. As expected (Sect. [J!2] i, the parameters of the surface- 
brightness profiles show larger errors for the bulge than for the 
disk. We recovered larger values for the Sersic parameters (r^., 
«) when the PSF FWHM is overestimated, and lower values 
when it is underestimated. Relative errors are correlated with 
the values of effective radius and shape parameter of the bulge 
but not with the magnitude of the galaxy. In Run 4 we built 
the artificial galaxies by adopting the correct PSF FWHM and 
a sky level that was 1 % larger than the actual one. For brighter 
galaxies an improper sky subtraction mostly affects the param- 
eters of the disk surface-brightness profile. For fainter galaxies, 
the large relative errors on the bulge parameters are due to their 
coupling with the disk parameters. This is consistent with the 
results of similar tests performed by Byun & Freeman (1995). 

The structural parameters to be measured to derive the in- 
trinsic shape of bulges are the ellipticity of the bulge and posi- 
tion angles of the bulge and disk (Sect. lS!2] i. In all the runs, the 
relative errors are smaller than 0.05 (absolute value) for galax- 
ies with Jt < 10 mag. Larger errors ( up to about 0.1) were 
found for fainter galaxies after an improper subtraction of the 
sky level from the image. However, this is not the case for our 
sample galaxies. 
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Table 1. Relative errors on the photometric parameters of the bulge and disk calculated for different galaxy magnitudes by means 
of Monte Carlo simulations. 









8 < < 9 






9 < Jt < 10 






10 < /t- < 1 1 






1 1 < Jj- < 12 




Parameter 


Run 


# 


Mean 


St.Dev. 


# 


Mean 


St.Dev. 


# 


Mean 


St.Dev. 


# 


Mean 


St.Dev. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


h 


1 


226 


-1 ■ 10"' 


0.02 


217 


-8 ■ I0-' 


0.04 


238 


-1 • 10-2 


0.06 


206 


-4 ■ 10-2 


0.14 




2 


223 


1 ■ 10^2 


0.08 


211 


1 ■ ur' 


0.09 


235 


9 • 10-2 


0.10 


205 


2 ■ 10-2 


0.16 




3 


225 


-1 ■ 10-2 


0.09 


214 


-I ■ ur' 


0.09 


231 


-8 ■ 10-2 


0.10 


208 


-6 ■ 10-2 


0.15 




4 


212 


-6 ■ 10-2 


0.08 


206 


-7 ■ 10-2 


O.IO 


220 


-6 ■ 10-2 


0.11 


180 


-3 ■ 10-2 


0.15 




1 


226 


1 ■ IQ-' 


0.02 


217 


7 ■ H)-' 


0.03 


238 


1 • 10-2 


0.06 


206 


6 ■ 10-2 


0.13 




2 


223 


-7 ■ 10-2 


0.08 


211 


-9 ■ 10-2 


0.08 


235 


-7 • 10-2 


0.09 


205 


-2 ■ 10-2 


0.14 




3 


225 


7 ■ 10-2 


0.06 


214 


9 ■ 10-2 


0.07 


231 


1 ■ 10-' 


0.09 


208 


1 ■ ur' 


0.13 




4 


212 


7-10-2 


0.09 


206 


8 ■ 10-2 


0.10 


220 


1 ■ 10-' 


0.11 


180 


8 ■ 10-2 


0.14 


n 


1 


226 


2 • 10-' 


0.02 


217 


4 ■ H)-' 


0.02 


238 


4 • 10-' 


0.04 


206 


4 ■ 10-2 


0.09 




2 


223 


-5 ■ 10-2 


0.05 


211 


-6 ■ 10-2 


0.05 


235 


-7 • 10-2 


0.07 


205 


-3 ■ 10-2 


0.12 




3 


225 


5 ■ 10-2 


0.04 


214 


6 ■ 10-2 


0.05 


231 


7 ■ 10-2 


0.07 


208 


1 ■ ur' 


0.11 




4 


212 


1 ■ 10-' 


0.08 


206 


1 ■ ur' 


0.08 


220 


1 ■ ur' 


0.09 


180 


1 ■ ur' 


0.15 


h 


1 


226 


-3 - 10--' 


0.03 


217 


-6 ■ I0-' 


0.03 


238 


-8 • 10-' 


0.05 


206 


-3 ■ 10-2 


0.14 




2 


223 


1 ■ 10-2 


0.05 


211 


1 ■ 10-2 


0.05 


235 


3 • 10-2 


0.07 


205 


2 ■ 10-' 


0.13 




3 


225 


-2 ■ 10-2 


0.05 


214 


-3 ■ 10-2 


0.05 


231 


-3 ■ 10-2 


0.05 


208 


-3 ■ 10-2 


0.07 




4 


212 


-7 ■ 10-2 


0.10 


206 


-9 ■ 10-2 


0.10 


220 


-9 • 10-2 


0.11 


180 


-1 ■ 10-' 


0.12 


h 


1 


226 


2 • 10-' 


0.02 


217 


9 ■ 10-' 


0.03 


238 


9 • 10-' 


0.04 


206 


4 ■ 10-2 


0.10 




2 


223 


3 • 10-' 


0.02 


211 


5 ■ 10-' 


0.03 


235 


6 • 10-' 


0.05 


205 


2 ■ 10-2 


0.11 




3 


225 


3 ■ 10-' 


0.02 


214 


8 ■ 10-' 


0.03 


231 


1 ■ 10-2 


0.05 


208 


6 ■ 10-2 


0.10 




4 


212 


8 ■ 10-2 


0.09 


206 


8 ■ 10-2 


0.08 


220 


9 ■ 10-2 


0.10 


180 


1 ■ ur' 


0.11 


Oh 
'io 


1 


226 


-1 ■ 10-' 


0.01 


217 


-4 ■ 10-' 


0.02 


238 


-2 • 10-' 


0.03 


206 


1 ■ 10-2 


0.06 




2 


223 


3 ■ 10-' 


0.01 


211 


3 ■ 10-' 


0.01 


235 


3 • U)-' 


0.02 


205 


1 ■ 10-2 


0.06 




3 


225 


-5 ■ 10-' 


0.01 


214 


-7 ■ 10-' 


0.02 


231 


-5 • 10-' 


0.03 


208 


-5 ■ 10-' 


0.06 




4 


212 


1 ■ 10-2 


0.10 


206 


4 ■ 10-2 


0.12 


220 


4 ■ 10-2 


0.12 


180 


6 ■ 10-2 


0.13 




1 


226 


-1 ■ 10-' 


0.01 


217 


-8 ■ 10-' 


0.04 


238 


-3 ■ 10-' 


0.03 


206 


-3 ■ 10-2 


0.09 




2 


223 


8 ■ 10-' 


0.02 


211 


8 ■ 10-' 


0.04 


235 


2 ■ 10-2 


0.04 


205 


2 ■ 10-2 


0.11 




3 


225 


-1 ■ 10-2 


0.02 


214 


-2 ■ 10-2 


0.03 


231 


-2 ■ 10-2 


0.04 


208 


-5 ■ 10-2 


0.10 




4 


212 


-8 ■ 10-2 


0.08 


206 


-7 ■ 10-2 


0.08 


220 


-9 • 10-2 


0.09 


180 


-1 ■ 10-' 


0.10 


PAb 


1 


226 


1 ■ 10-* 


0.01 


217 


1 ■ K)-' 


0.02 


238 


-2 ■ 10-' 


0.02 


206 


-4 ■ 10-' 


0.06 




2 


223 


-1 ■ 10-' 


0.02 


211 


-1 ■ ur' 


0.02 


235 


-2 ■ 10-' 


0.03 


205 


3 ■ 10-' 


0.05 




3 


225 


1 ■ 10-' 


0.01 


214 


3 ■ 10-' 


0.02 


231 


-3 ■ K)-' 


0.03 


208 


-I ■ K)-' 


0.06 




4 


212 


-1 ■ 10-2 


0.10 


206 


-2 ■ 10-2 


O.IO 


220 


-3 ■ 10-2 


0.10 


180 


-I ■ 10-2 


0.12 


PAd 


1 


226 


-3 ■ 10-' 


0.02 


217 


-3 ■ IO--* 


0.02 


238 


-2 ■ 10-' 


0.02 


206 


5 ■ 10-' 


0.07 




2 


223 


-1 ■ 10-' 


0.04 


211 


-1 ■ ur' 


0.02 


235 


1 ■ 10-' 


0.04 


205 


6 ■ 10-' 


0.08 




3 


225 


2- 10-' 


0.04 


214 


-2 ■ ur' 


0.04 


231 


-2 ■ 10-' 


0.04 


208 


5 ■ 10-' 


0.08 




4 


212 


8 ■ 10-' 


0.08 


206 


3 ■ 10-' 


0.08 


220 


4 ■ 10-' 


0.08 


180 


2 ■ 10-2 


0.11 



NOTE. Col.(l): Photometric parameter; Col. (2): Run of the Monte Carlo simulation. Artificial galaxies are analyzed by assuming the conect values of PSF FWHM and sky 
level (Run 1), con'ect sky level and a PSF FWHM 2% larger with respect to the actual one (Run 2), con'ect sky level and a PSF FWHM 2% smaller with respect to the actual 
one (Run 3). Artificial galaxies are built by assuming the conect PSF FWHM and a sky level 1% larger with respect to the actual one (Run 4); Cols. (3, 6, 9, 12): Number of 
artificial galaxies in the magnitude bin; Cols. (4, 7, 10, 13): Mean of the relative eiTors; Cols. (5, 8, 1 1, 14): Standard deviation of the relative eiTors. 



3.5. Results and comparison with previous studies 

The parameters derived for the structural components of the 
sample galaxies are collected in Table [3| All the listed val- 
ues are corrected for seeing smearing and galaxy inclina- 
tion. Surface brightnesses were calibrated by adopting, for the 
2MASS images, the flux zero point given in the image headers 
( i Jarre tt et al. 2000i) . 

The comparison of the structural parameters obtained for 
the same galaxy by different authors is often not straightfor- 
ward on account of possible differences in the observed band- 
pass, parameterization of the surface-brightness distribution, 
and fitting method. 

MHOl already studied 11 of our sample galaxies 
(NGC 772, NGC 2775, NGC 2841, NGC 2985, NGC 3169, 
NGC 3626, NGC 3675, NGC 3898, NGC 4450, NGC 4501 
and NGC 4826). They performed a two-dimensional paramet- 



ric decomposition of the /-band surface brightness distribu- 
tion. They considered ellipticities and position angles of both 
the bulge and disk as free parameters. Therefore, we considered 
their results as the most suitable to be compared with ours. The 
structural parameters we measured are consistent with those 
given by MHOl, within 25% for all the common galaxies but 
NGC 4826. We argue that they strongly underestimated the 
scale length of its disk (and consequently obtained a wrong 
estimate of the other parameters) because of the small field of 
view (3' X 3') of their image. In Fig.|2]we show the comparison 
between our axial ratios and position angles of the bulge and 
disk and those measured by MHOl. 

4. Correlations between structural parameters 

The study of correlations between the structural parameters of 
bulges and disks of our sample galaxies will help us to both 
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This is in agreement, within the errors, with the correla- 
tion by MHOl. If we use the mean surface brightness inside 
one effective radius instead of the effective surface brightness 
this relation becomes the so-called Korm endy relation, already 
known for bulges and elliptical galaxies ( iKormendv 1977 ). 

Finally, the absolute luminosity of the bulge is correlated 
with the effective radius. Larger bulges are more luminous (Fig. 
[3j-'). This result in 

Mb = -21.93(±0.06)- 3.4(+0.2) log re (r = -0.80), (11) 

where Mb is the 7-band magnitude of the bulge. A similar re- 
sult was obtained by MHOl for a smaller sample of disk galax- 
ies spanning a larger range of Hubble types. 

Bulges and elliptical galaxies follow a tight relation, the 
FP, defined by the effective radius, mean surface bright- 
ness within effective radiu s ^ and central veloci ty dispersion 



Fig. 2. Comparison between the axis ratios (left panel) and the 
position angles (right panel) measured in this paper and by 
MHOl. Filled dots and open diamonds correspond to values 
measured for bulges and disks, respectively. Residuals hq and 
APA are defined as 1 - qumxiqow and 1 - PAmhoi/PAqui-, re- 
spectively. 

cross check our results with those available in literature and 
identify and rule out peculiar bulges from any further analysis. 

4.1. Bulge parameters 

We did not find any correlations between the bulge parame- 
ters and Hubble type. Neither the effective radius (Fig. |3K), 
effective surface brightness (Fig. [3ji) nor the n shape parame- 
ter (Fig.[3p) show a statistically significant Pearson correlation 
coefficient (r). 

From near-infrared observations of spiral galaxies, 
Andredakis et al. (1995) found that bulges of early-type spirals 
are characterized by « a: 4 (i.e., they have a de Vaucouleurs 
radial surface brightness profile), while the bulges of late-type 
spirals are characterized by n » 1 (i.e., they have an expo- 
nential radial surface brightness profile). This early result 
was confirmed in various studies (e.g., de Jong et al. 1996; 
Khosroshahi et al. 2000; Mac Arthur et al. 2003; Mollenhoff 
& Heidt 2001; Mollenhoff' 2004; Hunt et al. 2004). We argue 
that our data does not show such a correlation due to the 
smaller range of Hubble types covered by our sample (SO-Sb) 
with respect to the cited works, where it is mostly evident for 
Hubble types later than Sb. 

The n shape parameter increases with effective radius. 
Larger bulges have a surface-brightness radial profile that 
which is more centrally peaked than that of the smaller bulges 
(Fig.[3p)). We obtained 



(IDjorgovski & Davis 19871 : iDressler et al. 19871) . Therefore, 



we derived the FP for the bulges of our sample galax- 
ies. The measurements of the central stellar velocity disper- 
sion for a subsample of 98 galaxies were available in litera- 
ture and were retri eved from the on line HyperLeda catalog 
(Patu rel et al. 20031) . The velocity dispersions given by the cat- 
alogue are corrected to a circular aperture of radius of 0.595 
h"' kpc, which is equivalent to an angular diameter of 3'.'4 at 
the distance of Coma, following the prescription by Jorgensen 
et al. (1995). The aperture-corrected velocity dispersions are 
given in Table [3] The coefficients describing the FP 

logre = 1.08(+0.09)logcro-i-0.25(+0.02)<;/e)-6.61(+0.40),(12) 

were derived by minimizing the square root of the residuals 
along the log axis. Errors given for every coefficient were 
calculated by performing a bootstrap analysis with 1000 it- 
erations. No statistically significant difference was observed 
when only bulges of lenticular or early-to-intermediate spiral 
galaxies were considered. The dispersion around this relation is 
cr = 0.1 1 dex and was measured as the rms scatter in the residu- 
als of log Tg. The observational error on the FP is 0.066 dex and 
includes the measurement errors in log (0.055 dex), log o"o 
(0.029 dex), and (//g) (0.021 mag arcsec"-). Errors in logrg and 
(jig) are not independent (Kormendy 1977). Compared with the 
dispersion around the relation, this gives an intrinsic scatter of 
0.088 dex. Figure |4] shows an edge-on view of the FP. Our co- 
efficients and those by Falcon-Barroso et al. (2002) are consis- 
tent within the errors, although they analyzed /T-band data. 
Unfortunately, we have not found in the literature the coef- 
ficient of the FP in the 7-band for a direct comparison (see 
Bernardi et al. 2003b . 

One of the projections of the FP is the so-called Faber- 
Jackson relation (FJ), which relates the luminosity of ellip- 
tical galaxies and bulg es to their central velocity dispersion 
(Faber & Jac kson 1976h . We derived the 7-band FJ relation for 
the bulge subsample obtaining 



log « = 0.38(+0.02)-H 0.1 8(+0.05) logre (r = 0.28). (9) log cro = 0.1(+0.2) - 0.095(±0.009)Mb (r = -0.71). (13) 



The effective surface brightness is dependent on the effec- 
tive radius. Larger bulges have a lower effective surface bright- 
ness (Fig.[3ji). We found a linear regression 

logA^e = 17.74(+0.07)H- 1.7(+0.2) logre (r = 0.55). (10) 



This result also holds when we consider only galaxies with 
errors on the central velocity dispersion lower than 10 km s~^ 
(A log (To < 0.018 dex). In fact, we derived 

log cro = 0.4(+0.3) - 0.085(+0.013)Mb (r = -0.65), (14) 
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Fig. 3. Correlations between the bulge parameters. Correlation between the Hubble type and effective radius (A), effective surface 
brightness (B), and shape parameter (C). Correlations between the effective radius and shape parameter (D), effective surface 
brightness (E), and absolute magnitude (F). In each panel the solid line represent the linear regression through all the points. The 
Pearson correlation coefficient (r), and the results of the Unear fit (y = a + bx) are also given. 



which is consistent within errors with Eq.[T3] From Eq.[T3]we 
derived L cc cr^^, which is very close to the virial relation and 
indicates that our bulges share important characteristics with 
bright elliptical galaxies dMatkovic & Guzman 2005h . On the 
other hand, Balcells et al. (2007) found L oc ctq^ (close to faint 
ellipticals) observing a sample of bulges with the Hubble Space 
Telescope in the /T-band. This discrepancy is not due to the 
adopted fitting method and it is not observed if only bright 
bulges (Mb < -20 mag) are considered. Indeed, we found 
L oc cr-^'^^^'^^ for our sample and L oc cr'^'^'^^'''^ for the sam- 
ple by Balcells et al. (2007). The different behaviour of faint 
bulges (Mb > -20 mag) requires further investigation to be 
explained. 

Khosroshahi et al. (2000) noticed that the shape parame- 
ter, effective radius and central surface brightness of elliptical 
galaxies and bulges are correlated. This relation was termed 



photometric plane (PP). Figure|6]shows an edge-on view of the 
PP of our bulge sample 

logn = 0.17(±0.02)logre-0.088(+0.004);uo+1.48(±0.05).(15) 

The coefficients were derived by minimizing the square 
root of the residuals along the log« axis. Errors given for 
every coefficient were calculated by performing a bootstrap 
analysis with 1000 iterations. No statistically significant dif- 
ference was observed when only bulges of lenticular or early- 
to-intermediate spiral galaxies were considered. The dispersion 
around this relation is cr = 0.04 and was measured as the rms 
scatter in the residuals of logn. Our coefficients and those by 
MHOl are consistent within the errors, although their sample is 
dominated by bulges of late-type spirals. The presence of the 
bulges of lenticular and spirals on the same PP hints towards a 
common formation scenario. 
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Fig. 4. Edge-on view of the FP for the 98 early-to-intermediate 
type bulges of our sample with measured velocity dispersion. 
The solid line represents the linear fit to the data. The dotted 
lines represent the Icr deviation from the fit. The eiTor bars in 
the lower right corner of the panel indicate the mean errors of 
the data. 



Fig. 5. FJ relation for the 98 early-to-intermediate type bulges 
of our sample with measured velocity dispersion. The solid line 
represents the linear fit to the data. The dotted lines represent 
the Icr deviation from the fit. The error bars in the lower right 
corner of the panel indicate the mean errors of the data. 



4.2. Disk parameters 

Regarding the disk parameters, we found no correlation be- 
tween the scale length and Hubble type (r = -0.06, Fig. IlK)- 
The same is true for the central surface brightness. In fact, it 
shows a large scatter also with Hubble type (r = -0.05, Fig. 
|7j}). This is consistent with the results of de Jong et al. (1996) 
and MHOl. 

On the other hand, the central surface brightness and the lu- 
minosity of the disks are dependent on the scale length. Larger 
disks have a lower central surface brightness (Fig. |7p). We 
found a linear regression 

log^io = 17.36(±0.1) -I- 1.4(+0.2) log/z (r = 0.49), (16) 
and brighter disks show larger scale lengths (Fig.|7][r) 

Md = -21.21(+0.09)-3.5(+0.2)log/! (r = -0.80), (17) 

where Md is the 7-band magnitude of the disk. 

The coefficients are in agreement within the errors with 
those given by MHOl. 

4.3. Bulge and disk interplay 

We have found that the disk scale length increases with central 
velocity dispersion. Since central velocity dispersion correlates 
with the virial mass of the bulge (Aih - ar^cr^jG with a - 
5, Cappellari et al. 2006), we conclude that larger disks are 
located in galaxies with more massive bulges (Fig.|8]\). For the 



subsample of 98 early-to-intermediate bulges with a measured 
velocity dispersion we found 

logcro = 2.13(±0.03)-h0.27(±0.06) \ogh (r = 0.42).(18) 

We also found a strong correlation between the bulge ef- 
fective radius and the disk scale length. Larger bulges reside 
in larger disks (Fig. [8}J). This relation was already observed 
by Courteau et al. (1996) and later observed in NIR by MHOl, 
MacArthur et al. (2003). We obtained a linear regression 

log re = -0.45(+0.03)-h0.91(±0.07) log/z (r = 0.74), (19) 

which is in agreement within error bars with the correlation 
found by MHOl. All these correlations between bulge and disk 
parameters indicate a link between the bulge and disk forma- 
tion and evolution history. This connection was interpreted as 
an indication of the for mation of late-type bu lges via secu- 
lar evolution of the disks ( Courteau et al. 19961) . However, our 
measurements of the scale lengths of the bulge and disk and 
n (Fig. [St) are also fully consistent with the predictions of 
the numerical simulations by Scannapieco & Tissera (2003) 
and Tissera et al. (2006). They studied the effects of merg- 
ers on the mass distribution of bulges and disks of galaxies 
formed in hierarchical clustering scenarios. Our mean value 
(re//i) = 0.36 + 0.17 is also in good agreement with rg/Zi found 
by Naab & Trujillo (2006) for a series of major mergers' rem- 
nants. These results indicate that these relations are not enough 
to distinguish between bulges formed by mergers or by secular 
evolution of the disk, even if a strong crosstalk between both 
components is present. 
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Fig. 7. Correlations between the disk parameters. Correlations between the Hubble type and disk scale length (A) and central 
surface brightness (B). Correlation between the disk scale length and absolute luminosity (C) and central surface brightness (D). 
Solid lines and coefficients as in Fig. [3] 



5. The equatorial intrinsic ellipticity of bulges 



5.1. Geometrical formalism 



In the previous section, we realized that the knowledge of cor- 
relations between structural parameters of the bulge and disk is 
not sufficient to distinguish between the different scenarios that 
were proposed to explain the formation of bulges. Therefore, 
we decided to study the intrinsic shape of bulges in order to 
give a further constraint to these scenarios. 

Independently of its internal structure, we can consider a 
bulge of a spiral galaxy as an ellipsoidal stellar system located 
in the center of the galaxy, which stands out against the disk in 
the photometric observations. We assume that both the bulge 
and disk share the same center, which coincides with the galac- 
tic one, and they have the same polar axis (i.e., the equatorial 
plane of disk coincides with that of bulge). 



Let {x,y,z) be Cartesian coordinates with the origin in the 
galaxy center, the jf-axis and y-axis corresponding to the prin- 
cipal axes of the bulge equatorial ellipse, and the z-axis corre- 
sponding to the polar axis. As the equatorial plane of the bulge 
coincides with the equatorial plane of the disk, the z-axis is 
also the polar axis of the disk. If A, B, and C are lengths of the 
ellipsoid semi-axes, the corresponding equation of the bulge in 
its own reference system is given by 



X 

A2 



y 

B2 



(20) 



Let {x',y',z') be now the Cartesian coordinates of the ob- 
server system. It has its origin in the galaxy center, the polar 
z'-axis along the line-of-sight (LOS) and pointing toward the 
galaxy, and the plane of the sky lies on the {x',y') plane. 
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Fig. 8. (A) Correlation between the disk scale length and central velocity dispersion for the 98 galaxies of our sample with a 
measured velocity dispersion. (B) Correlation between the disk scale length and bulge effective radius. (C) The ratio between 
the bulge and disk exponential scale lengths as a function of the bulge shape parameter. Filled circles and crosses represent the 
results of our measurements and simulations by Tissera et al. (2006), respectively. Solid lines and coefficients as in Fig.|3] 



The projection of the disk onto the sky plane is an ellipse 
whose major axis is the line of nodes (LON), i.e., the inter- 
section between the galactic and the sky planes. The angle 9 
between the z-axis and z'-axis corresponds to the inclination 
of the bulge ellipsoid; it can be derived as = arccos (cZ/c) from 
the length c and d of the two semi-axes of the projected ellipse 
of the disk. We defined (p (0 < <p < n/l) as the angle between 
the ;ic-axis and the LON on the equatorial plane of the bulge 
(x,y). Finally, we also defined ifr (0 < ij/ < n/T) as the angle 
between the x'-axis and the LON on the sky plane (x',y'). The 
three angles 0, </>, and t// are the usual Euler angles and relate 
the reference system (x, y, z) of the ellipsoid with that {x',y',z') 
of the observer by means of three rotations. Indeed, since the 
location of the LON is known, we can choose the x'-axis along 
it, and consequently it holds that i/' = 0. By applying these two 
rotations to Eq.|20]it is possible to derive the equation of the el- 



lipsoidal bulge in the reference system of the observer, as well 
as the equation o f the ellipse corr esponding to its projection 
on the sky plane dSimonneau et al . 1998). Now, if we identify 
this ellipse with the ellipse that forms the observed ellipsoidal 
bulge, we can determine the corresponding axes of symmetry 
Xp and y^.. The first one, on which we measured the semi-axis 
fl, forms an angle 6 with the LON (the x'-axis in the observed 
plane); the semi-axis b is taken over the y^-axis. We always 
choose < 6 < n/l, so it is possible that a either be the major 
or the minor semi-axis, and vice versa for b. 
We have 



a'h- = A^C^siifecos-<!> + B^Csitf0sia'ip + A-B^cos-g. 
~ + h~ = A"(cos" + cos" sin" 0) + S^(sin" + cos" (9cos" 1^) + C" sin" 
(B^ -A2)cosSsin2(^ 



tan 2S = 



A^fcos^Ssin i*-cos^i4) + B2(cos^ flcos^.*- sin ) + sin 



(21) 

(22) 

(23) 
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Fig. 6. Edge-on view of the PP for early-to-intermediate bulges 
of our sample. The solid line represents the linear fit to the data. 
The dotted lines represent the Icr deviation from the fit. The 
error bars in the lower right corner of the panel indicate the 
mean errors of the data. 



If the ellipsoidal bulge is triaxial {A + B C) then it is 
possible to observe a twist {6 + 0; see Eq.|23]l between the axes 
of the projected ellipses of the bulge and the disk. 

5.2. Inverse problem or deprojection 

We will now focus our attention on the inverse problem, i.e., 
deprojection. Following Simonneau et al. (1998), from Eqs.lSTl 
I22I and|23] we are able to express the length of the bulge semi- 
axes (i.e. A, B, and C) as a function of the length of the semi- 
axes (i.e. fl, b) of the projected ellipse and the position angle 
(5). 



a^+b'- 



1 + e cos26 + sinlS- 



, sine/) 1 

.1 

cos(p cosO 

SOS(j) 1 

7 

sin(^ cosO ^ 

. cosO cos2(f> 1 + cos^9 

1 + e I IsinlS — coslo 

sin'-Q sin2(f> sin'-O 



1 + e cos2S - sin26- 



(24) 
(25) 
(26) 



where e - {a^ - b^)/{a^ + b^) is, in some way, a measure of the 
eUipticity. It will be - 1 < e < I. 

Notice that a, b, 6 and 9 are all observed variables. 
Unfortunately, the relation between the intrinsic and projected 
variables also depends on the spatial position of the bulge (i.e., 
on the (p angle), which is not directly accessible to observa- 
tions. For this reason, only a statistical determination can be 
performed to assess the intrinsic shape of bulges. 

As A and B are the semi-axis of the equatorial ellipse of the 
bulge, we have to distinguish between two cases, according to 



Eqs.|24]and|25] If a > b (or equivalently e > 0) then A > B. 
Otherwise, if a < b (or equivalently e < 0) then A < B. Thus, if 
6 the equatorial plane of the bulge ellipsoid is not circular 
and the bulge ellipsoid is triaxial. 

From Eqs.|24l|25]and|26]we can write the axial ratios A/C 
and B/C as explicit functions of (p. Moreover, we assume that 
the angle (p is random and independent of the length of ellipsoid 
semi-axes. Thus, the normalized probability distribution P(0) 
of getting a given value of (p in ((p, (p+d(p) is 



Picp)dcp^\. 

Jo 



(27) 



According to a fundamental theorem of statistics, the prob- 
ability of obtaining a given value of any function f((p) (e.g., 
one of the axial ratios) will be equal to the probability of get- 
ting the corresponding value of <p, provided that the ratio df/d(p 
between the corresponding differential elements is taken into 
account. 

In this work, we will focus our attention on the intrin- 
sic equatorial eUipticity of bulges. We define it as £ = (A^- 
B^)/{A^+B^) with -!<£'< 1. In a forthcoming paper, we will 
study the intrinsic flattening of the bulge ellipsoids defined as 
the ratio between the length C of polar semi-axis and the mean 
length of the equatorial semi-axes. 

From Eqs.|24]and|25] it is straightforward to derive a rela- 
tion among the intrinsic variables (equatorial eUipticity E and 
position angle 0), and the measured ( i.e. 0, e, and 6), which is 



1 



e sin 26 



E sin(20) 
1 + £■ cos(20) cos 6* 1 + e cos 2(5 



= 2- 



(28) 



The second member of the equality in Eq. |28] allows us to 
define the observable Q in terms of the measured variables 6, 
e, and S. It must be stressed that for each specific bulge the 
relation between the equatorial eUipticity E and the unknown 
parameter (p embraces the whole of the measured variables 
through the single variable Q. 

On the one hand, Eq. |28] will yield the conditional proba- 
bility Pq(E) that a given bulge, with a measured value of Q, 
takes on any particular value of E (individual statistic); on the 
other hand, this equation will give the probability Pe{Q) asso- 
ciated to each value of Q for a bulge with intrinsic equatorial 
eUipticity E. This latter probability will be the kernel of an in- 
tegral equation that relates the observed statistical distribution 
P{Q), corresponding to a sample of galaxies, with the statisti- 
cal distribution of the equatorial eUipticity P(E) for the same 
sample. 

5.3. Individual statistics 

For a given galaxy, we can measure the values of 0, e and 6, 
and then derive the value of Q through Eq. |28] We want to de- 
termine the probability Pq(E) that such a galaxy (i.e. with such 
a value of Q) will take on a value of E in the range {E, E+dE). 
The subindex Q specifies this galaxy. All the galaxies with the 
same value of Q shall partake the same probability distribution 
Pq(E). 
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Fig. 9. Probability distribution function of E for the galaxy IC 
4310. The dotted line represents the value of T=0.098 derived 
for this galaxy. The value of £1/2 is also shown in the plot. 

Once the value of Q is prescribed, for some values of E 
there are not values of (0 < <^ < 7t/2) that satisfy Eq. |28] 
Hence, it shall hold that Pq(E) - 0. Only for those values of E 
such that 



E'> 



0<T <1, 



will there will exist two values of (p that fulfill Eq. |28] 

Then for any value of E > T the probability Pq{E) will be 
given by 



By calculating the partial derivative of Eq. [ 
izing we obtain 



PAE) = 



E sfgryp 
arccos T 



(30) 

]and normal- 
(31) 



where the subindex T plays the same role as the subindex Q\ 
both are related through Eq. |29l 

This means that the possible values of E are very concen- 
trated and slightly larger than T . To get an idea of how PjiE) 
is peaked near the value of T, we calculated the value £1/2 for 
which the total probability that E > £1/2 is equal to the prob- 
ability that E < Ei/2- For every bulge £1/2 is a sort of mean 
value of E, and is given by 



^1/2 



1 +T 



(32) 



In Fig|9]we show, as an example, the probability Pt{E) for 
one galaxy of our sample. 

5.4. Global statistics 

Likewise, we can define the probability Pe(Q) associated to 
each value of Q for a given bulge of intrinsic equatorial ellip- 
ticity E. 



For a prescribed value of E, it will only be possible to get 
the values of that satisfy Eq.|28]when T < E. The probability 
P{4>j) for the two values of (pj is given by Eq. |27] The proba- 
bility Pe(T) is equal to the sum of the two probabilities Pi<pj), 
weighted with the ratio {6(I>I6T)e of the differential elements: 



e(G) = X 



5Q 



Once the partial derivative from Eq. 
obtain 



Pe{T) = - 



1 



E>T 



(33) 



I is computed, we 



(34) 



Our purpose here is to determine the probability distribu- 
tion P(E), that is the number of galaxies with a value of E 
inside the interval (E, E + dE), starting from a sample that is 
sufliciently representative. We have measured for such a sam- 
ple the distribution P{T), namely the number of bulges whose 
values of T are within T and T + dT. We must write now P(T) 
as the integral over all the values of E of the product of the con- 
ditional probability Pe{T) by the so far unknown distribution 
P(E): 



(29) P(T) 



P(E)PE(T)dE. 



Then, by making use of Eq. 
tegral equation 



P(T) 



2 r' 

Jt 



P(E) 



V£2 - r2 



-.dE, 



(35) 

we obtain an Abel-like in- 
(36) 



which will allow us to derive P(E) from the observed distribu- 
tion P(T). 

However, as usual, the data PiT) of our statistical prob- 
lem takes the form of histograms, hence, the relevant equations 
must be formulated accordingly. 

Let Tk with = 0, 1,2,...,A' (To = 0, Jjv = 1) be a set 
of discrete ordinates that defines the histogram of the observed 
function PiT) (Fig. [TO]l. The A;'''-element of this histogram is 
defined by 



Pk(T) 



1 



Tk - Tk- 



r 



P(T) dT. 



(37) 



We must now seek the integral equation that relates the vari- 
ables Pk(T) with the probability distribution P{E). 

We notice that the integral of P{T) in Eq. [32]is equivalent 
to the difference between the two quadratures of PiT) over the 
intervals (Tk-i,!) and (Tk,l). In both of them, we will replace 
P{T) by its integral form, given by Eq.[36l 

Then, by inverting the order of integration, we can easily 
rewrite the two resulting integrals to obtain 



Pk(T) 



1 



Tk - Tk-\ 

2 In 
Tk - Tk-i 



P(E) dE 



(M(Tk-i)-M(Tk)), 



(38) 
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where we defined 



M(Tk) 



x: 



P{E) dxcsXniTklE) dE, 



(39) 



and M{Tk-\) in a similar way. 

Equation[38]is the integral equation that will allow us to ob- 
tain the values of P{E). Since it is consistent with the numerical 
structure of the data, we are confident that we have eliminated 
most of the numerical problems that arise from the direct inver- 
sion of Eq. |36] which constitutes a typical ill-posed problem. 

At this point, we require for P{E) the same histogram rep- 
resentation that we have already introduced for the data P{T). 
We introduce a similar set of discrete ordinates for the vari- 
able E (i.e. Ek ~ Tk), and an analogous definition to obtain the 
elements PkiE) of the histogram 



Pj(E) 



1 



r 



P{E) dE. 



Thus Eq.|38]can be rewritten into the form 



Pk{T) = PkiE) ■ 



Tk - Tk- 



■ {M(Tk-i) - M(Tk)) . 



(40) 



(41) 



In order to express the integrals M{Tk-i) and M(Tk) as lin- 
ear functions of the so far unknown values of Pj{E) (where 
7 > A; - 1), we consider that P{E) is constant and equal to the 
unknown values of Pj{E) over each interval (Ej^i, Ej), accord- 
ing to its histogram representation. Therefore, it is straightfor- 
ward to compute the coefficients CM{k, j), defined by the linear 
relation 



M(Tk)^ J] CMik,j)Pj(E). 



(42) 



=k+l 



Thus Eq. |4T] becomes a simple linear algebraical equation 
that relates the terms of the two histograms Pk{T) and Pj(E) 
through a triangular matrix. Once we have the integral Eq. [36] 
into a suitable matrix form, according to the histogram rep- 
resentation of the data and results, a simple matrix inversion 
could be, in principle, enough to obtain the resulting Pj{E). 
However, such a procedure may add to the intrinsic difficulties, 
due to the lack of precision typical of the observational data, 
which naturally arise in the matrix inversion process; a catas- 
trophic mixture when dealing with an inverse problem. 

5.5. Inversion methods for the integral equation 

We have considered two different approaches to tackle the nu- 
merical problem. The first one is suggested by the method that 
leads us to the integral Eq.|4T]for the set of discrete elements 
Pk(T). We notice that there are two different terms in its right- 
hand side. The first one is the identity operator The second one 
is the difference between two integrals of PiT), multiplied by 
a kernel that is quickly decreasing. We can consider the former 
as the leading term, and treat the latter as a corrective term in a 
iterative perturbation method. Consequently, we can write Eq. 
EI] as 

2/k 



Pk(E) = PkiT) + 



Tk - Tk-i 



mTk-i) - M{Tk)) , 



(43) 



where we can determine P^(£) at the /'''-iteration making use 
of the form of P^"'(£') from the (/ - l)'''-iteration to compute 
the correction term [M{Tk-\) - M(Tk)]. As an initial guess, we 
consider PkiE) equal to Pk{T) at a zero-order approximation. 
This iterative process is repeated until convergence is achieved. 

From the data P(T) of the histogram shown in Fig. [10] we 
have obtained a satisfactory solution P{E) with a few number 
(5-10) of iterations. We will discuss later the physical quality 
of this solution, but we must recognize here the stability of the 
method. Actually, we have always achieved the same solution 
with a few number of iterations, starting from any trial initial 
distribution (namely, any zero-order approximation for Pk(E)). 
Moreover, the greater advantage of solving the integral equa- 
tion by means of an iterative perturbation method is the possi- 
bility to recover, according to Eq.[4T] the approximate diagram 
P'i^(T) that corresponds to any iterative solution P'kiE), and this 
yields a double check on the evolution and quality of results. 

However, when dealing with this kind of inverse problem, 
it may often happen in the practice that the results are correct 
mathematically, but not from the physical standpoint. In view 
of this difficulty, and in spite of the excellent quality of the 
foregoing iterative method, we wished to develop an alternative 
method of inversion for the integral Eq. [36] in order to double 
check the results. 

The other way to numerically treat the integral Eq. [36] 
comes from the analytical inversion 



r p(E)^ f 

Jt' Jt' 



■P{T)dT'. 



(44) 



Once we have this analytical form for the required solu- 
tion, we rewrite it for the histogram representation of P{E), by 
defining PkiE) as Pk{T) in Eq. [40] Again the difference be- 
tween the two integrals of P(T) show up. Now, we impose the 
histogram model of this P(T) to derive analytically the matrix 
elements that relate any Pk(E) to all the elements PkiT). The 
corresponding solutions obtained with the two methods are the 
same, allowing for the small differences due to round-off errors 
of the two different numerical algorithms employed. 

Now that we are confident of the reliability of both methods 
of inversion of the integral equation, we can came back to the 
aforesaid difficulties. 

When applying either method to the observed distribution 
PiT), in form of a histogram with bins PkiT) as shown in Fig. 
[To] we may obtain non-physical results, these are negative val- 
ues for some bins PkiE). This occurs due to the fact that in the 
frame of the adopted histogram representation for PkiT) and 
PkiE) and the associated numerical algorithm chosen to repre- 
sent the matrix operator for the integral Eq. [36] the measured 
values PkiT) cannot be the integral transform of any physical 
distribution PiE). However, another set of values PkiT) that are 
slightly different from the original, and consequently compati- 
ble with the observations, might satisfy the above requirement; 
i.e., it can be the integral transform of some physical PiE) and 
its inverse transform will solve our problem. 

These considerations claim a statistical regularization pro- 
cess, which can be achieved by considering the histogram PiT) 
to be the statistical mean of 1000 histograms, all of them com- 
patible with the observations according to Poisson statistics. 
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Fig. 10. PDF of T. The probability is normalized over 10 bins, 
which are geometrically distributed to cover the interval (0, 1). 
The width of the first bin is 0.03 and the width ratio of two con- 
secutive bins is 1.25. Error bars correspond to Poisson statis- 
tics. 

For each one of the 1000 possible realizations for Pk{T), we 
have obtained the corresponding Pk(E) by means of the in- 
version of the integral equation. The non-physical histograms 
Pk(E), i.e., those with some negative bins, were rejected. From 
the physical solutions, we have obtained the mean histogram 
and the corresponding error bars, as shown in Fig.[TT] The sta- 
tistical regularization process also allows us to estimate errors 
due to the possible lack of statistics in the sample. 

5.6. The probability distribution function of intrinsic 
ellipticities 

In Fig.[TT]we present the PDF of the bulge intrinsic ellipticities 
P{E). It was obtained by applying the procedure described in 
Sect. 153] using the PDF P{T) shown in Fig.[TOl The T values 
for each galaxy were calculated by means of Eqs. l28]and|29l 
from the measured values of e, 6 and 6. 

The PDF is characterized by a significant decrease of prob- 
ability for E<Q.Q1 (or equivalently Z?/A<0.93), suggesting that 
the shape of bulge ellipsoids in their equatorial plane is most 
probably elliptical rather than circular Such a decrease is 
caused neither by the lack of statistics (because in the regu- 
larization method we took into account the Poisson noise) nor 
by the width of the bins (because we tried different bin widths). 

We have calculated the average E value weighted with the 
PDF through 



{E} 



ZiPDF(Ei)*Ei 
T.iPDF{Ed 



(45) 



obtaining a value of {E) = 0.16 «B/A) = 0.85). This is fully 
consistent with previous findings by Bertola et al. (1991) and 
Fathi & Peletier (2003), based on the analysis of smaller sam- 
ples of bulges. For the sake of comparison, the value of (B/A) 
was derived from their data using Eq. 45. It is (B/A) = 0.85 
for the bulges studied by Bertola et al. (1991). They adopted a 
different approach to derive the PDF for intrinsic axial ratio of 



Fig. 11. PDF of E. The probability is normalized over 10 bins, 
which are geometrically distributed to cover the interval (0, 1). 
The width of the first bin is 0.03 and the width ratio of two 
consecutive bins is 1.25. The error bar of each Pk{E) bin cor- 
responds to the Poisson statistics of 1000 realizations of Pk(T) 
after excluding non-physical cases. 



bulges from the misalignment of the major axes of bulges and 
disks and the apparent ellipticity of bulges. It is (B/A) - 0.79 
for the early-type bulges of the sample studied by Fathi & 
Peletier (2003). They measured the bulge equatorial ellipticity 
by analyzing the deprojected ellipticity of the ellipses fitting 
the galaxy isophotes within the bulge radius. 

A further important result derived from P(E) is that there 
are not bulges with E > 0.6 (B/A < 0.5). This is also in 
good agreement with Bertola et al. (1991) and Fathi & Peletier 
(2003). They found a minimum axial ratio B/A = 0.55 and 
B/A = 0.45, respectively. 

We also studied the possible differences in the shape of 
bulges depending on their observational characteristics (i.e., 
morphology, light concentration, and luminosity, see Fig. fTST i. 
First of all, we subdivided our bulges according to the mor- 
phological classification -3 < HT < and < HT < 3 of 
their host galaxies to look for differences between lenticular 
and early-type spiral galaxies. A second test was done by sub- 
dividing the sample bulge between those with a Sersic index 
n < 2 and those with n > 2 to investigate possible correlations 
of bulge shape with light concentration. Finally, we subdivided 
our bulge into faint (Mb > -22) and bright (Mb < -22) in or- 
der to search for differences of bulge shape with the 7-band 
total luminosity. We did not find any significant difference be- 
tween the studied subsamples. They are characterized by the 
same distribution of E, as confirmed at a high confidence level 
(> 99%) by a Kolmogorov-Smirnov test. 

Several authors discussed the problem of the intrinsic shape 
of elliptical galaxies by means of observations and/or numeri- 
cal simulations. Ryden (1992), Lambas et al. (1992), and Bak 
& Statler (2000) agree that the observed distribution of elliptic- 
ities cannot be reproduced by any distribution of either prolate 
or oblate spheroidal systems. Any acceptable distribution of tri- 
axial systems is dominated by nearly-oblate spheroidal rather 
than nearly -prolate spheroidal systems. The formation of triax- 
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ial elliptical galaxies via simulation of merging events in the 
framework of a hierarchical clustering assembly was studied 
by Barnes & Hernquist (1996), Naab & Burkert (2003) and 
Gonzalez-Garcia & Balcells (2005). On the other hand, in the 
monolithic scenario where the galaxy formation occurs at high 
redshift after a rapid collapse, we may expect that the final 
galaxy shape would be nearly spherical or axisymmetric, as 
recently found in numerical experiments by Merlin & Chiosi 
(2006). But there is no extensive testing of the predictions of 
numerical simulations against the derived distribution of bulge 
intrinsic ellipticities. 

5. 7. The influence of nuclear bars on P{E) 

The presence of nuclear bars in galaxy bulges has been known 
since the former work of de Vaucouleurs (1974). However, it is 
only in the last decade with the advent of high-resolution imag- 
ing that a large number of them have been detected allowing the 
study of their demography and properties (see Erwin 2004, and 
reference therein). 

The sample galaxies were selected to not host large-scale 
bars, according to visual inspection and photometric decompo- 
sition of their 7-band images (Sect.|2]i. These selection criteria 
did not account for the presence of nuclear bars. In fact, our 
sample has 23 galaxies in common with the samples studied by 
Mulchaey & Regan (1997), Jungwiert et al. (1997), Martini & 
Pogge (1999), Marquez et al. (1999) and Laine et al. (2002). 
They were interested in the demography of nuclear bars. A nu- 
clear bar was found in 6 to 8 out of these 23 galaxies (26%- 
35%), according to the different authors' classifications. 

Since nuclear bars are more elongated than their host bulges 
and have random orientations, they could affect the measure- 
ment of the structural parameters of bulges and consequently 
their PiE). To address this issue we carried out a series of 
simulations on a large set of artificial galaxies. They were ob- 
tained by adding a nuclear bar to the artificial image of a typical 
galaxy of the sample and analyzing the structural properties of 
the resulting image with GASP2D, as done in Sect. | 3. II 

We adopted a Ferrers profile (ILaurikainen et al. 2005h to 
describe the surface brightness of the nuclear bar component 

2\«nb+0.5 



/nb(^,/7) 




if r„b < flnb 
if r„b > flnb 



where the coordinates (^, rf) are defined as in Sect. 13. II and flnb, 
/o,nb, and «nb are the bar length, its central surface brightness, 
and a shape parameter describing the curvature of the surface- 
brightness profile, respectively. Following Laurikainen et al. 
(2005) we chose Wnb = 2. Nuclear bar isophotes are ellipses 
centered on (^o, ^o) with constant position angle PAnb and con- 
stant ellipticity tnb = 1 - ?nb- The radius rnb is given by 

'•nb = [(-(^ -^o) sin PAnb + (?7-?7o) cos PAnb)^+ 

- ((^ - ^o) cos PAnb + {r]- m) sin ^Kbflqlb\'^ ■ (46) 

We generated 1000 images of galaxies with a Sersic bulge, 
an exponential disk, and a Ferrers nuclear bar. The structural 



parameters of the bulge and disk were selected to match those 
of a typical galaxy of the studied sample. It has Jj - 9.6 mag, 
re = 0.87 kpc, n = 2.32, h = 2.47 kpc and BIT = 0.37 ac- 
cording to the mean values of the structural parameters given 
in Table [3] Apparent ellipticity and position angle of the bulge 
and disk were randomly chosen in the ranges defined in Sect. 
13.41 The structural parameters of the nuclear bar were ran- 
domly chosen in the ranges < flnb < '"e for the length, 
0.2 < ^nb < 0.7 for the ellipticity, 0° < PA„b < 180° for the po- 
sition angle, and < L„b/T < 0.02 for the nuclear bar-to-total 
luminosity ratio. 

The flnb range was estimated from the 5 sample galaxies 
with a nuclear bar in common with Laine et al. (2002), which 
are characterized by (flnb/'"e) - 0.8. Detailed studies about 
luminosities of nuclear bars are still missing. Nevertheless, 
the LnhIT range was derived by considering that some nu- 
clear bars are secondary bars, which reside in large-scale bars. 
According to Erwin & Sparke (2002), a typical secondary bar 
is about 12% of the size of its primary bar. From Wozniak 
et al. (1995) we derived that the luminosity of the secondary 
bar is about 18% of that of the primary one. Since a primary 
bar contribu tes abo ut 15% to the tota l luminosity of its galaxy 
jPrieto et al. 20011: ILaurikainen et al. 2005h . the typical L^^IT 
ratio for a nuclear bar is about 2%. 

All simulated galaxies were assumed at a distance of 30 
Mpc. Pixel scale, CCD gain and RON, seeing, background 
level and photon noise of the artificial images were assumed 
as is Sect. 13.41 The two-dimensional parametric decomposition 
was applied by analyzing with GASP2D the images of the ar- 
tificial galaxies as if they were real. We defined errors on the 
parameters as the difference between the input and output val- 
ues. The mean errors on the fitted ellipticity and position angle 
of the bulge and disk and their standard deviation are given in 
Table |2] They correspond to systematic and typical errors. 

Table 2. Errors on the ellipticity and position angles of the 
bulge and disk calculated for galaxies with nuclear bars by 
means of Monte Carlo simulations. 



Parameter 


Mean 


St.Dev. 




0.02 


0.03 




0.005 


0.02 


APAb n 


0.3 


7 


APAd n 


0.1 


3 



For each sample galaxy the values of e, S, and 9 were de- 
rived in Sect. |5] from the ellipticity and position angle mea- 
sured for the bulge and disk. We randomly generated a series 
of q<o, qd, PAb, and PAd by assuming they were normally dis- 
tributed with the mean and standard deviation given in Table |2] 
We tested whether bulges are axisymmetric structures, which 
appear elongated and twisted with respect to the disk compo- 
nent due to the presence of a nuclear bar. To this aim, for each 
galaxy we selected 1000 realizations of q^,, q^, PAb, and PAj 
which gave smaller e (i.e., a rounder bulge) and smaller S (i.e., 
a smaller misalignment between bulge and disk) with respect 
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Fig. 12. PDF of E for the different subsamples. (A) Lenticular galaxies (-3 < HT < 0). (B) Spiral galaxies (0 < HT < 3). (C) 
Galaxies with n > 2. (D) Galaxies with n < 2. (E) Faint bulges (Mb > -22). (F) Bright bulges (Mb < -22). Bin widths and error 
bars as in Fig. [TT] The mean intrinsic ellipticity and the number of galaxies of the subsample are given in each panel. 



to the observed ones. This correction can be considered as an 
upper limit to the bulge axisymmetry. We obtained 1000 PiT) 
distributions and calculated P(E) from their mean. 

Following this procedure, if we consider that all galaxies 
in our sample host a nuclear bar, we obtain a P{E) where the 
decrease of the probability for E < 0.07 disappears (Fig. [T3T l. 
This means that most of the bulges are circular in the equato- 
rial plane. The average value of the ellipticity of (£) - 0.12 
i{B/A) = 0.89). However, if we consider a more reahstic frac- 
tion of galaxies that host a nuclear bar (i.e, 30% as found by 
Laine et al. 2002 and Erwin & Sparke 2002), the resulting P{E) 
is consistent within errors with the P(E) derived in Sect. 15.61 
(Fig.inil. The average value is {E} = 0.15 ({B/A} = 0.86). 

6. Conclusions 

The structural parameters of the bulge and disk of a magnitude- 
limited sample of 148 unbarred SO-Sb galaxies were inves- 
tigated to constrain the dominant mechanism at the epoch of 
bulge assembly. 

- We presented a new fitting algorithm (GASP2D) to perform 
two-dimensional photometric decomposition of galaxy im- 



8 



6 - 




0.0 0.2 0.4 0.6 0.8 1.0 

E 

Fig. 13. PDF of E for the original sample (solid lines), for a 
sample with 30% of bulges with a nuclear bar (dotted line) and 
for a 100% fraction of galaxies hosting a nuclear bar (dashed 
line). Bin widths and error bars as in Fig.fTTI 

ages. The surface-brightness distribution of the galaxy was 
assumed to be the sum of the contribution of a Sersic bulge 
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and an exponential disk. The two components were char- 
acterized by elliptical and concentric isophotes with con- 
stant (but possibly different) elhpticity and position angles. 
GASP2D is optimized to deal with large image samples, 
and it adopts a robust Levenberg-Marquard fitting algo- 
rithm in order to obtain reliable estimates of the galaxy 
structural parameters. 

- The bulge and disk parameters of the sample galaxies were 
derived from the y-band images, which were available in 
the Two Micron All Sky Survey. 

- The bulges of the sample galaxies follow the same FP, FJ, 
and PP relationships found for elhptical galaxies. No statis- 
tically significant difference is observed when only bulges 
of lenticular and early-to-intermediate spiral galaxies were 
considered. This supports the idea that bulges and elhpti- 
cals formed in the same way. 

- Tight correlations between the parameters of bulges and 
disks were found. In fact, the disk scale lengths increase 
with both the central velocity dispersion and bulge effec- 
tive radius. Therefore, larger disks reside in galaxies with 
more massive and larger bulges. This was interpreted as an 
indication of the formation of bulges via secular evolution 
of their host disks. 

- Our measurements of the exponential scale length of the 
bulge and disk, as well as of bulge shape parameter, were 
also fully consistent with numerical simulations of the ef- 
fects of mergers on the mass distribution of the bulge and 
disk in galaxies formed in hierarchical clustering scenarios. 

- These results indicate that the above relations are not 
enough to clearly distinguish between bulges formed by 
early dissipative collapse, merging or secular evolution. 
All these mechanisms could be tested against the intrinsic 
shape of bulges. Therefore, the PDF of the intrinsic equa- 
torial elhpticity of the bulges was derived from the distri- 
bution of the observed elhpticities of bulges and their mis- 
alignments with disks. 

- About 80% of bulges in unbarred lenticular and early-to- 
intermediate spiral galaxies are not oblate but triaxial el- 
hpsoids. Their mean axial ratio in the equatorial plane is 
(B/A) = 0.85. This is consistent with previous findings by 
Bertola et al. (1991) and Fathi & Peletier (2003). There 
is no significant dependence of the PDF on the morphol- 
ogy, Ught concentration, and luminosity of bulges. The de- 
rived PDF is independent of the possible presence of nu- 
clear bars. 
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2.43 


0.84 


1.5 


18.06 


18.30 


0.87 


6.7 


NGC5292 


(R')SA(rs)ab 


9.62 


4729 


63.0 


- 


17.09 


3.6 


1.51 


0.65 


51.3 


17.82 


15.1 


0.81 


49.3 


NGC5313 


Sb 


9.80 


2732 


37.8* 


- 


16.15 


1.5 


0.62 


0.72 


33.2 


17.42 


12.6 


0.51 


45.8 


NGC5326 


SAa: 


9.80 


2712 


37.8* 


165 


18.17 


9.1 


1.65 


0.53 


133.4 


18.93 


22.3 


0.60 


133.0 


NGC5440 


Sa 


9.83 


3890 


51.9 


- 


17.53 


3.4 


1.95 


0.70 


41.7 


18.42 


15.7 


0.36 


46.1 


NGC5485 


SAO pec 


9.32 


2074 


32.8* 


159 


18.18 


8.9 


2.27 


0.85 


0.9 


17.93 


15.3 


0.78 


173.4 


NGC5533 


SA(rs)ab 


9.75 


4051 


54.0 


- 


17.70 


5.1 


2.34 


0.66 


23.6 


18.03 


14.1 


0.64 


26.3 


NGC5614 


SA(r)ab pec 


9.50 


4073 


54.3 


- 


17.60 


6.1 


3.41 


0.94 


47.2 


17.84 


12.8 


0.88 


149.3 


NGC5631 


SA(s)0° 


9.39 


2075 


32.7* 


168 


17.29 


6.1 


3.61 


0.82 


125.4 


17.55 


11.6 


0.94 


157.7 


NGC5687 


SO" 


9.89 


2192 


34.4* 


190 


18.25 


7.5 


2.31 


0.66 


103.7 


19.72 


27.1 


0.54 


96.9 


NGC5838 


SAO- 


8.48 


1555 


28.5* 


266 


17.38 


6.2 


2.01 


0.82 


44.9 


18.65 


30.6 


0.34 


42.8 


NGC6278 


SO 


9.96 


2786 


37.1 


150 


16.83 


3.1 


2.35 


0.66 


109.7 


18.08 


12.6 


0.53 


131.4 


NGC6340 


SA(s)0/a 


9.25 


1169 


22.0* 


144 


17.17 


4.3 


3.17 


0.95 


121.5 


17.59 


16.0 


0.97 


64.6 


NGC6851 


SO 


9.69 


2913 


38.8 


224 


17.31 


5.9 


3.13 


0.69 


167.4 


17.81 


10.9 


0.73 


156.7 


NGC6861 


SA(s)0-: 


8.66 


2636 


35.5* 




16.67 


9.1 


2.11 


0.53 


142.4 


18.01 


17.8 


0.70 


133.8 



Table 3: continued. 
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Jt 




D 
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n 
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h 


Id 
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(km s"') 


(Mpc) 


(kms-i) 
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(°) 
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o 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


NGC6890 


SA(rs)ab 


9.99 


2323 


31.8* 


- 


18.69 


5.4 


3.80 


0.55 


24.8 


17.18 


9.8 


0.70 


162.3 


NGC6893 


SA(s)° 


8.95 


2996 


39.9 


- 


16.86 


5.6 


1.85 


0.72 


2.1 


18.74 


26.4 


0.57 


10.9 


NGC6920 


SA(rs)0°: 


9.39 


2762 


36.8 


- 


16.54 


3.9 


2.46 


0.82 


140.4 


17.59 


13.1 


0.84 


130.2 


NGC7007 


SAO": 


9.91 


2801 


37.3* 


125 


16.78 


3.1 


1.86 


0.81 


133.0 


17.80 


11.7 


0.60 


1.7 


NGC7020 


(R)SA(r)0+ 


9.77 


2987 


37.8* 


- 


17.52 


2.8 


2.56 


0.89 


160.6 


18.20 


15.3 


0.40 


163.9 


NGC7049 


SA(s)0 


8.14 


2051 


27.6* 


240 


16.92 


8.4 


2.17 


0.74 


59.9 


17.04 


18.1 


0.76 


55.8 


NGC7083 


SA(s)bc 


9.37 


2999 


38.7* 


71 


18.19 


3.6 


1.34 


0.95 


16.1 


17.88 


18.6 


0.56 


7.9 


NGC7096 


SA(s)a 


9.84 


2847 


36.7* 


248 


18.95 


12.5 


4.05 


0.90 


124.8 


18.58 


10.8 


0.82 


101.3 


NGC7135 


SAO- pec 


9.72 


1799 


34.7* 


- 


18.84 


8.1 


4.69 


0.74 


6.6 


18.90 


21.1 


0.64 


48.8 


NGC7166 


SAO- 


9.47 


2233 


30.2* 


- 


16.12 


2.7 


1.70 


0.67 


9.9 


17.73 


13.5 


0.47 


13.2 


NGC7172 


Sapec 


9.44 


2315 


33.9* 


179 


17.35 


3.3 


1.16 


0.82 


80.3 


17.56 


15.5 


0.52 


95.8 


NGC7192 


SAO": 


9.37 


2794 


31.6* 


179 


18.08 


10.0 


2.91 


0.96 


8.5 


18.38 


16.5 


0.96 


137.0 


NGC7302 


SA(s)0-: 


9.99 


2221 


33.7* 


- 


16.31 


2.2 


1.35 


0.86 


88.0 


17.66 


11.0 


0.61 


99.0 


NGC7311 


Sab 


9.96 


4185 


55.8 


- 


16.11 


1.8 


1.26 


0.82 


25.8 


17.14 


9.3 


0.51 


12.2 


NGC7377 


SA(s)0+ 


9.08 


3036 


40.5 


144 


17.93 


7.7 


1.78 


0.84 


105.6 


18.22 


21.6 


0.78 


102.8 


NGC7550 


SAO- 


9.93 


4789 


63.8 


234 


18.68 


9.2 


3.40 


0.97 


179.5 


18.30 


12.1 


0.86 


146.8 


NGC7585 


(R')SA(s)0-"pec 


9.36 


3104 


41.4 


219 


17.74 


6.9 


3.33 


0.69 


106.4 


17.70 


15.1 


0.82 


106.4 


NGC7600 


SO- sp 


9.82 


3086 


41.1 


210 


17.96 


5.9 


2.75 


0.97 


175.0 


17.89 


11.7 


0.90 


150.1 


NGC7606 


SA(s)b 


8.62 


1881 


28.9* 


147 


17.45 


2.8 


1.02 


0.63 


147.2 


18.22 


31.9 


0.41 


145.1 


NGC7625 


SA(rs)a pec 


9.90 


1272 


23.0* 


- 


17.77 


4.3 


0.89 


0.59 


58.1 


16.91 


9.0 


0.74 


32.0 


NGC7702 


(R)SA(r)0+ 


9.92 


3062 


40.8 




16.64 


2.6 


1.02 


0.61 


99.8 


17.82 


12.9 


0.53 


114.8 


NGC7711 


SO 


9.99 


3705 


49.4 


180 


18.62 


5.7 


3.30 


0.75 


93.0 


19.37 


19.1 


0.33 


9 1.7 


NGC7722 


so/a 


9.88 


3673 


49.0 


165 


17.26 


4.3 


1.52 


0.61 


144.8 


17.37 


10.5 


0.85 


139.8 


NGC7742 


SA(r)b 


9.60 


1294 


22.2* 


95 


15.66 


1.4 


3.10 


0.91 


12.6 


16.24 


7.5 


0.98 


180.0 


NGC7769 


(R)SA(rs)b 


9.90 


3877 


51.6 




18.43 


6.6 


3.82 


0.60 


11.6 


17.59 


12.1 


0.87 


100.4 


UGC12591 


so/a 


9.99 


6615 


88.2 


288 


16.78 


3.7 


1.51 


0.43 


59.2 


17.77 


11.3 


0.51 


58.2 



NOTE. Col. (1): Galaxy name; Col. (2): Morphological classification from de Vaucouleurs et al. (1991, hereafter RC3); Col. (3): Total observed J-band magnitude from the 2MASS 
all-sky extended source catalogue (XSC); Col. (4): Radial velocity with respect to the CMB radiation from Lyon Extragalactic Database (Hereafter LEDA); Col. (5): Distances obtained 
as Vi/i/H() with Ho=75 km s-'Mpc-' , distances marked with an asterisk have been taken from Tully (1988); Col. (6): Central velocity dispersion from LEDA; Col. (7): Effective surface 
brightness of the bulge; Col. (8): Effective radius of the bulge; Col. (9): Shape parameter of the bulge; Col. (10): Axis ratio of the bulge; Col. (11): Position angle of the bulgescale 
length of the disk; Col. (12): Central surface brightness of the disk; Col. (13): Scale length of the disk; Col. (14): Axis ratio of the disk; Col. (15): Position angle of the disk. 



